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摘要 : 林 火 预测 预报 是 科学 有 效 进行 林 火 管理 的 前 提 , 是 林业 管理 部 门 和 科研 工作 者 的 育 泛 关 深 的 领域 。 人 逻辑 斯 带 回归 
(Logistic Regression, LR) 是 目前 国内 外 广泛 应 用 于 森林 火灾 预测 的 模型 方法 ,然而 近年 张 有 学 者 发 现 该 方法 没有 充分 考虑 林 
火影 响 因子 的 空间 相关 性 和 异 质 性 ,从 而 导致 模型 拟 合 结果 偏差 。 地 理 加 要 逻辑 斯 蒂 回 归 ( Geographically weighted logistic 
regression, GWR) 模型 考虑 到 了 模型 变量 之 间 的 空间 相关 性 ,有 效 提 高 的 模型 的 拟 合 能 力 。 为 探讨 GWLR 模型 在 福建 林 火 预 
测 上 的 适用 性 ,本 研究 应 用 LR 和 GWLR 两 种 方法 分 别 建立 福建 省 森林 火灾 与 气象 因子 的 预测 模型 ,通过 模型 拟 合 能 力 对 比 ， 
判断 在 GWLR 的 适用 性 。 研 究 以 2000 一 2005 年 福建 地 区 森林 火灾 江 遇 火 点 数据 和 每 日 气象 因子 为 基础 ,将 全 样本 分 为 60% 的 
建 模 数据 和 40% 的 校 验 数据 ,并 重复 5 次 ,建立 5 个 样本 组 。 选 择 在 5 个 样本 组 中 3 个 及 以 上 表现 显著 的 变量 进入 最 终 模型 。 


变量 进 
研究 结果 表明 GWLR 在 模型 拟 合 度 、 模 型 残 差 空间 自 相关 性 以 及 预测 准确 率 等 方面 均 优 于 LR 模型 ,说 明 充 分 考虑 模 


2 


型 变量 


的 空间 异 质 性 有 助 于 提高 模型 的 预测 精度 ,同时 也 验证 了 GWLR 在 福建 地 区 林 火 预测 上 的 适应 性 。 此 外 ,模型 参数 结果 显示 ， 
“日 最 高 地 表 气温 ” “日 最 低地 表 气 温 ” “日 平均 风速 ”\ 2 丽人 小 时 降水 量 ” “日 最 高 本 站 气压 “日 照 时 数 ”“ 日 最 高 气温 ”和 


“日 最 小 相对 湿度 ”8 个 因子 对 福建 省 林 火 发 生 有 显著 影响 ,研究 结论 为 福建 地 
关键 词 : 林 火 预测 ;空间 异 质 性 ;逻辑 斯 带 回 归 虽 t 理 加 权 逻 辑 斯 带 回 归 


区 林 火 预测 预报 提供 了 新 的 方法 。 
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Abstract: Forest fire" forecasting is a key component of effective and science-based forest management and has been 
cómprehensively addressed in the scientific literature. The logistic regression ( LR) technique has been used in forest fire 
prediction models. However, some scholars have recently reported that the technique does not adequately consider the 
spatial correlation and heterogeneity of fire impact factors , which results in poorly fitting models. In contrast, geographically 
weighted logistic regression ( GWR) models consider the spatial correlation of model variables, which improves the model's 


goodness of fit. In order to explore the applicability of the GWLR model in Fujian forest fire forecasting, the present study 
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used both the LR and GWLR methods to establish forecast model for forest fires and meteorological factors in Fujian 
Province, and the model fitting ability of two models were compared. Based on the forest fire and meteorological data for 
Fujian from 2000 to 2005, the original dataset was randomly divided into training (6096) and validation (4096) samples, 
with five replications and five sample groups, and predictors that were significant ( a — 0.05) for at least three of the five 
sample groups were included in the final models. The goodness of fit, residual error, spatial autocorrelation, and prediction 
accuracy of the GWLR model were all better than those of the LR model, and the GWLR comprehensively explained, the 
spatial heterogeneity of model variables and helped to improve the prediction accuracy of the model. The study also. verified 
the suitability of the GWLR model on the forest fire forecasting in Fujian area. In addition, the results also indicated that the 
occurrence of Fujian forest fires is significantly affected by eight parameters, including minimum and maximum surface 
temperature, daily average wind speed, daily precipitation, highest station pressure, hours of sunshine, daily maximum 
temperature, and daily minimum relative humidity. Therefore, the GWLR model may provide a/new technique for the 


prediction of forest fires in Fujian Province. 


Key Words: forest fire forecast; spatial heterogeneity; logistic regression; geographically, weighted logistic regression 
( GWLR) 


M JC e BRUM: ds 0 I RES s] DD ^ , X RP EUR E RB Xa bs 1. 2 RAK HE: i M 0 CARIUS, 
判定 林 火 发 生 的 影响 因子 ,建立 准确 的 林 火 预测 预报 模型 对 林 火 管理 卫 作 至 关 重 要 。 目 前 ,国内 外 有 关 林 火 
发 生 的 预测 模型 和 影响 因子 分 析 已 大 量 开 展 , 主要 模型 方法 为 人 工 神经 网 络 '” ERRES SERE, 
泊 松 回归 、 负 二 项 回归 FERAH . 零 膨 胀 负 二 项 以 及 梁 辑 斯 幕 回归 模型 "等 ,其 中 逻辑 斯 带 模 型 应 用 最 
为 广泛 。 然 而 以 上 所 有 模型 均 假 设 因 变量 与 自 变 量 之 间 的 关系 是 空间 平稳 的 , 即 模 型 参数 在 整个 研究 区 域 上 
是 一 个 不 变 的 常数 ,通常 建 模 的 结果 是 一 套 模型 参数 应 用 于 整个 研究 区 域 ,也 称 之 为 全 局 模型 。 然 而 , 随 着 研 
究 的 深入 ,很 多 学 者 发 现 区 域 的 空间 异 质 性 是 不 答 忽 略 的 。 林 火 与 影响 因子 之 间 的 关系 在 空间 上 是 非 平稳 
的 ,具有 很 强 的 异 质 性 所 ,因此 ,以 往 的 全 局 模型 在 林 火 预测 预报 上 可 能 会 产生 较 大 误差 。 

地 理 加 权 回 归 模 型 (GCWR) 考 虑 和 插 地 理 空间 因素 的 影响 ,该 模型 将 一 个 大 的 数据 集 分 成 了 若干 小 区 域 , 减 
少 了 各 个 模型 间 的 差异 性 ,有 助 于 提高 模型 精度 ,可 以 用 来 解决 空间 的 平稳 性 问题 。 目 前 ,地 理 加 权 回 归 模 型 
主要 应 用 在 生态 学 .经 济 学 和 社会 科学 等 领域 [5 ”1 。 近 年 来 ,国外 已 有 少数 学 者 将 该 模型 的 拓展 模型 即 地 理 
加 权 逻 辑 斯 蒂 模 型 (CWLR9 咬 用 于 森林 火灾 的 预测 预报 及 林 火 影响 因子 的 空间 分 析 上 ,并 且 表 现 出 较 好 的 拟 
合 效果 ”中 。 由 于 该 模型 强调 月 变量 与 因 变 量 关系 的 空间 异 质 性 ,因此 模型 在 不 同 区 域 的 适用 性 可 能 有 较 
大 差异 。 为 了 探讨 GWLR 模型 对 我 国 亚 热带 地 区 林 火 预测 预报 的 适用 性 , 且 由 于 2000 一 2005 年 福建 省 的 林 
火 发 生 既 不 是 太 多 ,也 不 是 很 少 ,比较 有 代表 性 ,因此 本 文 以 福建 省 2000—2005 年 卫星 火 点 数据 为 基础 ,选取 
每 日 气象 数据 为 月 变量 因子 ,对 模型 的 拟 合 能 力 进行 分 析 ,并 与 传统 全 局 逻辑 斯 蒂 回 归 模 型 进行 对 比 ,探讨 地 
理 加 权 回 归 模 型 在 林 火 预测 预报 上 的 适用 性 。 


1 材料 和 方法 


1. 了 0 研究 区 域 概况 

福建 省 位 于 我 国 东 南 沿海 地 区 , 介 于 北纬 23"33' 一 28"20' 东经 115°50' 一 120°40' 之 间 ( 图 1) 。 根 据 第 八 
次 全 国 森 林 资 源 清 查 结果 显示 ,福建 省 林业 用 地 面积 926.82 万 hm ,和 森林 面积 801.27 万 hm ,人 工 林 面积 
377.69 万 hm? ,森林 覆盖 率 65.9596 , 森林 蓄积 量 60796.15 万 m ,是 我 国 南方 重点 林 区 的 省 份 之 一 。 但 是 ， 
也 是 我 国 森林 火灾 的 高 发 区 。1951 一 1987 年 ,全 省 共 发 生 6.2 万 次 的 森林 火灾 ,受害 山林 面积 有 150 万 hm’, 
平均 每 年 森林 火灾 0.1 万 次 ,受害 山林 面积 4 万 hm? , 占 全 省 年 平均 造林 面积 的 2596— 33.3967" ;1998 一 2007 
年 间 ,福建 省 共 发 生 森 林 火 灾 4504 起 ,火场 面积 累计 达 79,572hm ,受害 森林 面积 累计 达 55485.6 hm? , 共 损 失 
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T 97 万 mm 的 林木 蓄积 量 ,其 中 2004 年 发 生 了 1164 次 
火灾 ,1998 年 林 火 发 生 次 数 最 少 ,为 156 次 ; 仅 10 年 就 
RET 40 次 重大 火灾 '”i。 福建 省 森林 火灾 问题 尤为 
严峻 。 
1.2. 数据 来 源 与 处 理 
1.2.1 数据 来 源 

福建 省 2000—2005 年 的 林 火 发 生 数据 (包括 起 火 
时 间 ,起 火 原因 .起火 地 理 坐 标 等 ) 来 源 于 林业 科学 数 
JE 中 ò ( http://www. cfsde. org/indexAction. action? 
classId = 1) 提供 的 卫星 火 点 解 译 数 据 ; 历史 气象 数据 来 
源 于 中 国 气 象 数据 共 享 网 络 (http ://cdc.cma.gov.cn/) , 


60 120km 
= ny 


3g 2000—2005 年 福建 省 内 22 个 国家 级 气象 站 的 每 日 116° — 117 118 — 119 — 120 — I2I*E 
气象 数据 ,一 共 包 22 个 气象 因子 。 由 于 仪器 设备 和 气 图 1 研究 区 域 示 意图 
象 站 点 变化 等 原因 ,造成 部 分 数据 缺失 , 且 无 法 修补 。 Po de 


KE , AIEO] CAR ER EE TT TOL TE , AILES S e EC 
据 , 剩 下 的 气象 因子 包括 日 平均 地 表 气 温 (% ) ,日 最 高 地 表 气 温 (% ) , Bois I REG (00) ,日 平均 风速 (my 
s) ,日 最 大 风速 (m/s) ,24 小 时 降水 量 ( mm) ,日 平均 本 站 气压 (hPa) ; 目 最 高 本 站 气压 (hPa) ,日 最 低 本 站 气 
压 (hPa) ,日 照 时 数 (hour) ,日 平均 气温 (人 ) ,日 最 高 气温 (CO 日 最 低 气 温 (%C ) ,日 平均 相对 湿度 (%) ,日 最 
小 相对 湿度 (% ) 共 15 个 气象 因子 。 
1.2.2 ”数据 处 理 

本 文 以 2000 一 2005 年 的 卫星 火 点 数据 为 基础 ,从 中 随机 选取 60% 的 训练 样本 (5210 个 林 火 数据 ) 来 构建 
模型 , 剩 下 40% 的 测试 样本 (3473 个 林 火 数据 为 用 来 检验 异型。 同时 ,为 了 避免 样本 分 布 对 试验 结果 的 影响 ， 
本 文 重复 5 次 试验 , 即 重复 $ 次 训练 样本 和 测试 样本 的 随机 划分 ,选择 5 次 试验 中 出 现 3 次 或 以 上 的 显著 变 
量 进入 全 样本 数据 拟 合 。 

本 文 分 别 应 用 SPSS 19.0 和 GWRA 软件 对 逻辑 斯 蒂 回 归 模 型 和 地 理 加 权 逻 辑 斯 蒂 回归 模型 进行 拟 合 , 并 
应 用 Rookcase 软件 对 模型 残 差 进行 空间 自 相 关 分 析 。 
1.8 ”研究 方法 
1.3.1 二 项 逻辑 斯 带 回 归 模 型 

(1) 多 重 共 线 性 QMulticollinearity ) 是 指 在 线性 回归 模型 中 , 自 变量 之 间 存 在 精确 的 相关 关系 或 者 高 度 相 
关 关 系 的 现象 。 

本 文 运用 方差 膨胀 因子 (the variance inflation factor , VIF) 诊断 法 对 林 火 发 生 的 驱动 因子 进行 共 线 性 检验 。 
方差 膨胀 因子 最 早 是 由 Marquardt 于 1960 年 引入 的 ,第 i 个 变量 的 方差 膨胀 系数 的 表达 式 为 


VIF, = 


lem (i - 1,2,...,n) (3) 


式 中 ,RR 为 自 变 量 x, 对 其 余 (n-1) 自 变 量 做 线性 回归 分 析 的 复 相 关系 数 。VIF, 的 值 越 大 , R? 越 接近 于 1 ,说 明 
第 i 个 自 变 量 与 其 他 (n-1) 个 自 变量 之 间 存 在 越 强 的 共 线 性 现象 ”1。 当 0<VIF<10 时 ,说 明 自 变量 之 间 
不 存在 多 重 共 线 性 现象 ; 当 10<VIF<100 时 ,说 明 自 变量 之 间 存 在 较 强 的 多 重 共 线性 现象 ; 当 VIF>100 时 ,说 
明 自 变量 之 间 存 在 严重 的 多 重 共 线性 现象 。 

(2) 惧 辑 斯 带 回归 模型 属于 广义 线性 回归 模型 ,对 于 只 有 两 种 分 类 结果 的 分 类 因 变 量 ,可 以 应 用 二 元 俱 
辑 斯 蒂 回 归 模 型 进行 建 模 分 析 , 其 预测 值 为 在 [0,1] 之 间 的 预测 概率 。 逻 辑 斯 蒂 模 型 是 目前 国内 外 应 用 最 广 
的 用 于 预测 林 火 发 生 概 率 的 模型 。 设 有 林 火 发 生 (Y= 1) 的 概率 为 已, 则 无 林 火 发 生 (Y= 0) 的 概率 为 (1- 忆 ) ， 
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则 有 林 火 发 生 的 概率 
e 1 
P(Y = 1) = 一 = 1 
( ) l+e l+e (1) 


WP , z 为 解释 变量 x ,x, ,… ,%, 的 线性 函数 ， 
z 20, OX, + Ox, t... + Ox, (2) 
式 中 ,ao ,ai ,0 ,... A, 为 各 个 解释 变量 的 回归 系数 。 
对 PP 进行 Logit 变换 ,即将 PP 变换 为 In[ P/(1-P) ], 则 有 
P 
]-P 
最 后 ,运用 极 大 似 然 估计 法 可 求 得 模型 的 参数 估计 系数 ao ,al ,as ,ay 77, 
1.3.2. 地理 加 权 逻 辑 斯 带 回归 模型 
由 于 地 理 位 置 的 不 同 而 引起 变量 之 间 的 关系 或 者 结构 发 生变 化 的 现象 称 为 空间 非 平 稳 性 (spatial 
nonstationarity) 。 地 理 加 权 逻 辑 斯 带 回归 模型 是 传统 逻辑 斯 带 回 归 模 型 的 扩展 , 示 态 了 空间 位 置 因素 ,利用 加 
权 最 小 二 乘法 对 每 一 个 坐标 点 进行 参数 估计 ,对 参数 的 估计 是 局 部 而 非 全 局 的 ,每 个 位 置 均 有 相应 的 参数 
估计 系 数 '*”]。 地 理 加 权 逻 辑 斯 带 回 归 模 型 (GWLR ) 的 表达 式 为 /与 逻辑 斯 蒂 模 型 相同 ,位 置 i 有 林 火 发 生 
(了 =1) 的 概率 为 P, 则 无 林 火 发 生 (Y=0) 的 概率 为 (1-P) , 则 位 置 i 有 林 火 发 生 概 率 


l 


=z =Q 十 QIX + aX, +... ta,x, (3) 


P(Y=1) exp( Ao) * 01,54 * Qo Xa Hla + 008,5) e 
1+ exp( a, 3 Q (us) 十 Oo i Xi He F O (ua) Yin ) l +e 
1 
EET. 
式 中 ,4— AO u; v) + O4 (v) Yil * AX vj) Yiz Te 十 Q hui vi) Vim o 
经 Logit 变换 ,有 

. P 

logit( P) = In "e. =ou 0) * Uis F Causa T + Oto) is 


最 后 ,运用 加 权 最 小 二 乘法 相 求 得 位 置 ; 的 局 部 回归 模型 的 参数 估计 系数 , 即 
&u;,v;) = (X W(u, v) X) "X'W(u, ,v;)logit( P) 

式 中 (u,,v;) 为 位 置 i BHO EAR BS, Y, s) 为 位 置 i 的 因 变 量 , 7 为 误差 项 。 Qu; ,v;) 为 模型 参数 Alu; vi) ， 
E EEA WS 的 估计 值 , W(u,;,v;) DEEE, X HMT EEE, EP , ERRE W 的 做 
法 为 :对 于 给 定 的 位 置 记 以 ;为 中 心 画 一 个 半径 为 > 的 圆 ,然后 计算 圆 内 每 个 观测 点 7 到 圆心 ; 的 距离 d; ,从 
而 生成 一 个 空间 权重 第 阵 W (00209992. AS paz H] GWLR 模型 对 火 点 进行 参数 估计 ,并 结合 ArcGIS 10.2 的 局 
部 多 项 式 插 值 法 对 整个 研究 区 域 进行 系数 插值 估计 1 。 
1.3.3 模型 预测 准确 率 计算 

本 文 根 据 ROC 曲线 分 析 法 所 计算 出 来 的 敏感 性 值 和 特异 性 值 , 可 求 得 约 登 指数 , 即 约 登 指数 = 敏感 性 值 
+ 特异 性 值 -1, 进 而 计算 出 福建 省 林 火 发 生 的 最 佳 临界 值 (cut-off 值 ) ,如 果林 火 发 生 的 预测 概率 大 于 该 临界 
值 则 认为 有 林 火 发 生 , 小 于 该 临界 值 则 认为 无 林 火 发 生 , 从 而 根据 模型 所 计算 出 来 的 林 火 发 生 次 数 与 实际 值 
进行 比较 ,进而 计算 出 LR 模型 和 GWLR 模型 模型 对 林 火 发 生 与 否 的 正确 判别 率 。 
1.3.4 ”模型 空间 自 相关 检验 Moran's I 

本 文 应 用 全 局 Moran's 工 指数 计算 残 差 ( 残 差 = 观测 值 -预测 值 ) 的 空间 自 相 关 来 对 模型 进行 评价 ,全 局 
Moran's 工 指数 的 值 越 小 ,说 明 残 差 的 空间 依赖 程度 越 低 ,模型 考虑 了 更 多 的 空间 结构 问题 ,模型 的 效果 越 好 。 

全 局 Moran's [指数 的 计算 公式 为 : 
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b (e; = e)? 
式 中 ,为 研究 区 域内 空间 单元 的 个 数 ,e 和 和 e; 分 别 为 空间 单元 ; 和 7 上 残 差 值 , e 为 残 差 e 的 平均 值 ,$ 为 残 差 
的 标准 差 ,o 是 空间 权重 矩阵 WW 中 对 应 的 元 素 ,如 果 空 间 单 元 i 和 j 相 邻 , 则 w = 1; 若 不 相 邻 , 则 w = 2 
Moran's I 指数 的 取 值 范围 为 [ -1,1] ,Moran's 工 指数 的 值 大 于 0, 说 明 研 究 区 域 的 残 差 值 呈 空间 正 相 关 ,小 于 0 
表示 残 差 值 呈 空间 负 相 关 , 等 于 0 则 表明 残 差 不 存在 空间 自 相关 性 2 。 


2 结果 与 分 析 


2.1 模型 拟 合 结果 与 分 析 
2.1.1 多 重 共 线 性 诊断 结 

多 重 共 线 性 诊断 结果 显示 ,变量 “日 平均 地 表 气 温 ”、“ 日 平均 本 站 气压 ”“ 日 最 低 本 站 气压 ”“ 日 平均 气 
温 ”“ 日 最 低 气 温 ” 等 5 个 变量 均 存 在 共 线 性 关系 ,将 这 5 个 变量 剔除 之 后 ,用 剩 下 的 Jo 个 自 变 量 构建 LR 模 
型 和 GWLR 模型 。 
2.1.2 LR 模型 拟 合 结 

本 文 应 用 LR 模型 和 “wald 向 前 ”原则 对 5 个 训练 样本 进行 拟 合 讳 算 , 得 到 5 个 不 同 的 特征 变量 子 集 , 然 
后 在 5 个 特征 变量 子 集中 选择 出 现 3 次 或 3 次 以 上 的 特征 变量 进 兴 全 样本 数据 的 拟 合计 算 ( 表 1) 。 


R1 LR 模型 特征 变量 选择 结果 


Table 1 LR model characteristic variable selection results 


变量 样本 1 aky — RURSÓ 样本 4 ”样本 5 最终 样 本 样本 显著 个 数 
Variable Sample 1 Sample 2 Sample 3 Sample 4 Sample 5 Final sample No.samples signif 
"E E 
日 最 再 地 表 气 温 十 十 十 十 十 十 3 
The highest surface temperature 
日 最 低地 表 气 温 
t 十 十 十 十 十 3 
The minimum surface temperature 
日 平均 风速 Daily average wind speed 十 十 十 十 十 十 5 
日 最 大 风速 Daily maxium wind speed 一 一 一 一 一 0 
24 小 时 降水 量 Daily precipitation 十 十 十 十 十 十 5 
最 高 本 站 气压 
日 最 再 本 站 气压 十 十 十 十 = + 4 
The highest station pressure 
日 照 时 数 Hours of sunshine 十 十 十 十 + + 5 
日 最 高 气温 Daily maximum temperature 十 十 十 十 十 十 5 
日 平均 相对 湿 刘 . E | 
Daily average'relative humidity 
日 最 小 相对 湿度 
十 十 十 十 ex 十 4 


The minimum relative humidity 


+ 表示 变量 在 模型 里 面 ;- 表 示 变 量 不 在 模型 里 面 


由 表 1 可 知 ， 日 最 高 地 表 气 温 ” “日 最 低地 表 气 温 "“ 日 平均 风速 “日 降水 量 ” “日 最 高 本 站 气压 ”、 
“日 照 时 数 ” “日 最 高 气温 ”和 ”日 最 小 相对 湿度 ”等 8 个 变量 进入 了 全 样本 数据 的 拟 合 阶段 , 且 除 了 ”日 最 高 
本 站 气压 ”和 “日 最 小 相对 湿度 ”在 5 个 中 间 模 型 特征 变量 选择 结果 中 出 现 4 次 以 外 ,其 他 6 个 变量 均 出 现 了 
5 次 。 表 2 为 8 个 显著 变量 的 全 样本 数据 拟 合 结果 ,结果 显示 各 变量 与 林 火 发 生 具 有 显著 相关 性 。 其 中 “日 
最 高 地 表 气 温 ” “日 平均 风速 ” “日 最 高 本 站 气压 “日 照 时 数 ” 和 ”日 最 高 气温 ”等 5 个 气象 因子 与 森林 火 
灾 呈 正 相 关 关系 “日 最 低地 表 气 温 ”“ 日 降水 量 ”" 和 ”日 最 小 相对 湿度 ”3 个 气象 因子 与 森林 火灾 呈 负 相关 关 
系 , 且 这 8 个 气象 因子 对 林 火 发 生 均 有 显著 的 影响 。 
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表 2 全 样本 数据 的 LR 模型 参数 拟 合 结果 

Table 2 LR model parameter fitting results of final sample 
变量 系数 标准 误差 显著 性 
Variable Coefficient standard error Significant 
日 最 高 地 表 气 温 The highest surface temperature 0.005084 0.000512 <0.0001 
日 最 低地 表 气 温 The minimum surface temperature —0.019823 0.000968 «0.0001 
日 平均 风速 Daily average wind speed 0.011462 0.001706 «0.0001 
24 小 时 降水 量 Daily precipitation —0.004977 0.000959 «0.0001 
日 最 高 本 站 气压 The highest station pressure 0.000257 0.000068 0.000157 
日 照 时 数 Hours of sunshine 0.012671 0.001224 <0.0001 
日 最 高 气温 Daily maximum temperature 0.006379 0.001082 «0.0001 
日 最 小 相对 湿度 The minimum relative humidity —0.007462 0.002545 0.003365 
常量 Constant -4.113521 0.701953 «0.0001 


2.1.3 GWLR 模型 拟 合 结 


首先 假设 模型 的 数据 结构 是 空间 非 平 稳 的 ,对 5 个 训练 样本 进行 模型 拟 合 ,然后 对 因 变 量 与 解释 变量 之 
间 的 关系 的 空间 非 平稳 性 进行 检验 。 如 果 “ 某 解释 变量 的 估计 系数 的 四 分 位 数 范 围 大 玫 LR 模型 中 该 解释 变 


量 的 估计 系数 的 + 1 标准 差 范 围 ”"”” , 则 认为 该 解释 
和 GWLR 模型 系数 估计 结果 见 附 表 1) 。 最 后 在 5 个 训练 样本 中 选择 出 现 3 次 或 3 次 以 上 的 非 平稳 变 


变量 是 显 


车 空间 非 平稳 变量 (以 样本 1 为 例 ,LR 模型 
入 


量 进 


全 样本 数据 的 拟 合 阶 段 , 并 作为 全 样本 数据 模型 的 非 平稳 项 ,而 其 他 平稳 变量 则 作为 全 样本 数据 模型 的 平稳 


结果 显 


项 。 检验 


示 除 了 “日 最 低 相对 湿度 "在 5 个 中 间 模 型 利 有 1 次 被 检验 为 空间 平稳 变 
象 因子 在 5 个 中 间 模型 中 均 被 检验 为 空间 非 平稳 变 


量 之 外 ,其 他 9 个 气 
量 ( 附 表 2)。 全 样本 数据 的 GWLR 模型 参数 拟 合 结果 如 


表 3 所 示 。 
表 3 全 样本 数据 的 GWLR 模型 参数 估计 
Table 3 Final sample of GWLR model parameter estimation 

系数 最 小 值 ”四 分 之 一 分 位 数 。 平均 数 中 位 数 ”四 分 之 三 分 位 数 最 大 值 
Coefficient Min 2596 quartile Mean Median 75% quartile Max 
ERIE Intercept 241.9643 -13.2123 -9.5616 -6.8187 -2.1540 2.8362 

E er ET 
d e ua —0.0023 0.0029 0.0040 0.0036 0.0050 0.0103 
The highest surface temperature 

zi DT 
H 最 低 电表 气温 一 0.0340 一 0.0264 一 0.0232 一 0.0234 一 0.0200 -0.0146 
The minimum surface temperature 
日 平均 风速 Daily average wind speed 一 0.0434 一 0.0036 0.0091 0.0100 0.0235 0.0591 
日 最 大 风速 Daily maxium wind speed —0.0261 —0.0106 0.0004 0.0008 0.0118 0.0274 
24 小 时 降水 量 

v 降水 量 -0.0173 -0.0071 -0.0054 -0.0050 -0.0031 0.0009 
Daily precipitation 
Hg f^ Ut 

RR US —0.0003 0.0001 0.0008 0.0006 0.0010 0.0040 
The highest station pressure 
日 照 时 数 Hours of sunshine —0.0005 0.0081 0.0112 0.0120 0.0145 0.0214 
日 最 高 气温 Daily maximum temperature 一 0.0014 0.0068 0.0119 0.0110 0.0159 0.0309 

平均 相对 温度 
ds 均 相 对 湿度 Ey —0.0883 -0.0217 -0.0092 -0.0016 0.0092 0.0244 
Daily average relative humidity 
日 最 小 相对 湿度 

; 相 湿度 D —0.0263 —0.0095 —0.0006 -0.0021 0.0055 0.0496 
The minimum relative humidity 

注 : 系 数 的 最 大 和 最 小 值 若 符号 相同 表示 该 变量 在 整个 研究 区 域 上 与 林 火 的 相关 性 一 致 ,反之 表示 变量 在 整个 区 域 上 与 林 火 的 相关 性 具有 

正 负 差 异 
表 3 显示 , 除 “ 日 最 低地 表 气 温 ” 在 整个 研究 区 域 上 均 与 林 火 发 生 呈 人 负 相 关 关 系 之 外 ,其 他 9 个 气象 因子 


在 整个 研究 区 域 上 均 在 正 相 关 和 人 负 相 关 之 间 变化 。 为 更 好 体现 GWLR 模型 各 变量 系数 的 局 部 变化 , 本文 运 
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用 ArcGIS 10.2 对 各 变量 模型 估计 系数 进行 空间 插值 (图 2)。 图 2 表明 GWLR 模型 的 估计 系数 随 空间 位 置 变 
化 而 变化 ,具有 明显 的 空间 异 质 性 。 


日 最 低地 表 气 温 
ISI -0.034—-0.030 
E -0.029—-0.026 


E -0.025—-0.022 
ES -0.021—-0.019 
IS -0.018—-0.014 


日 最 大 风速 
ES -0.026—-0.015 
E -0.014—-0.005 
E] -0.004 一 0.006 
E 0.007 一 0.016 
Wl 0.017 一 0.027 


日 最 高 本 站 气压 
IBS -0.0003—0.0005 
E 0.0006—0.0013 
E 0.0014—0.0022 
E 0.0023—0.0031 
I 0.0032—0.0040 


日 最 高 地 表 气 温 
IS -0.002 一 0 
E 0.001 一 0.002 
E] 0.003 一 0.005 
看 0.006 一 0.007 
I 0.008—90.102 


»z 


日 平均 风速 
IBS -0.043 一 -0.022 
E -0.021—0 

E 0.001—0.018 
IS 0.019—0.039 
Ill 0.040—0.059 


24 小 时 累计 降水 
E -0.017—-0.013 
E -0.012—-0.010 


|..170.009—-0.006 
IS -0.005—-0.003 
I -0.002—-0.001 


日 照 时 数 

mS -0.001 一 0.004 
E 0.005—0.008 
|..10.009—0.013 
E 0.014—90.017 
Bi 0.018—0.022 
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日 最 高 气温 
E -0.002—0.004 


日 平均 相对 湿度 
IR -0.087—-0.065 
E -0.064—-0.043 
[—]-0.042—-0.021 
E -0.020—-0.001 
IB 0.002—0.023 


常量 
EN -41.158—-32.371 
[73 -32.370—-23.585 


日 最 小 相对 湿度 


[E] 0.005—0.019 E -23 584—-14.799 
E 0.020—0.034 E -14.798—0 
EE 0.035—0.049 I8 0.001—2.774 


图 2 变量 系数 分 布 图 


Fig.2 Variable coefficient distribution 


IEAk, Speed Rt RR cS EGZET T8 [8] E, ARA ;检验 值 的 绝对 值 小 于 1.96 , 则 表示 估计 
系数 在 研究 区 域 上 不 显示 , 若 t HRAT- 1.96 或 大 于 1.96 ,说 明 估计 系数 显著 。 图 3 表明 模型 变量 的 估计 
系数 的 显著 性 也 具有 很 强 的 空间 异 质 性 、 

图 2, 图 3 综合 显示 “日 最 低地 表 气 温 ” 在 整个 研究 区 域 上 均 与 林 火 发 生 呈 负 相 关 关 系 , 且 该 因子 的 估计 
系数 在 整个 研究 区 域 P354: 24 ^ NES ERZK t" 与 林 火 发 生 呈 显著 负 相 关系 数 关 日 照 时 数 " 和 “日 最 高 气 
温 ” 与 林 火 发 生 旦 显著 正 相关 关系 ;“ 日 最 高 本 站 气压 ”与 林 火 发 生存 在 正 负 两 种 相关 性 ,但 以 正 相 关 为 主 , 主 
要 分 布 在 福建 北部 和 南部 地 区 ;日 平均 相对 湿度 "在 南平 市 北部 地 区 和 龙岩 市 大 部 分 地 区 与 林 火 呈 显 著 的 
负 相 关 关 系 , 仪 在 三 明 与 南平 市 交界 处 存在 小 块 正 相关 区 域 ;“ 日 最 小 相对 湿度 ”在 龙岩 市 西南 地 区 存在 显著 
正 相 关系 痊 i 在 三 明 市 和 永安 市 交界 处 及 厦门 市 和 漳州 市 的 小 部 分 沿海 地 区 存在 显著 的 负 相 关系 数 ;日 平 
均 风 速 ?得 福建 省 西部 和 中 部 地 区 与 林 火 发 生 呈 显著 正 相 关 , 而 “日 最 大 风速 "在 西部 和 中 部 地 区 与 林 火 发 生 
WEREWERE 
2.2 模型 评价 

本 文 将 全 样本 分 成 60% 的 训练 样本 (用 于 建 模 ) 和 40% 的 测试 样本 (模型 检验 )。 应 用 最 小 信息 准则 
C AIC) RŽI ASSE) F ROC 曲线 (ROC) 下 的 面积 AUC 值 (AUC ) 等 统计 量 和 模型 的 预测 准确 率 对 LR 
和 GWLR 模型 进行 对 比 评价 。 

模型 拟 合 统计 结果 表明 ,对 比 传统 逻辑 斯 蒂 回 归 模 型 ,地理 加 权 逻 辑 斯 蒂 回 归 模 型 具有 更 小 的 AIC 和 
SSE 值 ,和 更 大 的 AUC 值 以 及 更 高 的 模型 预测 准确 率 ( 表 4)。 表 4 显示 ,LR 模型 的 预测 准确 率 为 74.2% 一 
76.2% 小 于 GWLR 模型 (77.1% 一 78.6% ) 。 
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N 
A 
日 最 高 地 表 气 温 
- 高 : 6.791 
低 : 1.960 
日 最 低地 表 气 温 日 平均 风速 
高 : -5.002 高 : 5.461 
低 : -12.793 低 : 1.961 
高 : -1.960 
低 : -3.068 
日 最 大 风速 24 小 时 累计 降水 
高 : 3.952 高 : -1.960 
低 : 1.960 低 : -3.531 
高 : -1.960 
低 : -4.576 
日 最 高 本 站 气压 日 照 时 数 
高 : 6.245 m 高 : 7.272 
低 : 1.960 低 : 1.960 


E 高 : -1.962 
低 : -2.075 
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日 最 高 气温 日 平均 相对 湿度 
高 : 7.598 高 : 2.263 
" 低 : 1.960 低 : 1.960 
- 高 : -1.960 
低 : -4.983 
T^ wb o l4 
日 最 小 相对 湿度 m^ HE" 常量 
高 : 2.984 Wd xo 高 : -1.960 
低 : 1.961 低 : -6.608 
-= 高 : -1.961 
低 : -4.983 
0 60 120km 
图 3 变量 系数 显著 性 分 布 图 
Fig.3 Significant distribution of variable coefficient 
正 显著 相关 性 用 暖色 调 表示 , 负 显 著 相 关 性 用 冷色 调 表示 
X4 模型 拟 合 统计 量 
Table 4 Model fit statistics 
men! 残 关 平 方 和 wia 预测 正确 率 Prediction accuracy/96 
样本 模型 Akaike Sum of Squares 。 下 的 面积 临界 值 训练 样本 测试 样本 
Sample Model information for Error Area under Cut-off a Erud E 
criterion ( AIC) (SSE) Curve ( AUC) ee (4095) 
样本 1 LR 5256.368 872.679 0.822 0.501030 75.5 74.1 
Sample 1 GWLR 4999.887 785.656 0.857 0.524587 78.6 76.8 
样本 2 LR 5307.347 883.882 0.822 0.501080 75.2 75.2 
Sample 2 GWLR 5054.525 795.851 0.855 0.466608 78.5 76.9 
样本 3 LR 5304.786 883.694 0.819 0.479230 75.1 75.9 
Sample.3 GWLR 5116.377 822.763 0.843 0.513953 71.6 71.3 
样本 4 LR 5342.380 891.813 0.815 0.503125 74.9 75.5 
Sample 4 GWLR 5072.533 800.321 0.853 0.540291 77.8 771 
样 杰 5 LR 5401.440 905.102 0.810 0.501975 74.2 76.2 
Sample 5 GWLR 5172.060 820.554 0.845 0.522827 71.6 79.0 
全 样本 LR 8792.525 1464.812 0.821 0.502310 75.1 
Complete sample GWLR 8323.496 1328.274 0.853 0.505205 78.4 
2.3 IRIN 


根据 两 个 模型 对 5 样本 和 全 样本 数据 的 拟 合 结果 ,分别 绘制 模型 的 残 差 图 (图 4) 。 图 4 表明 GWLR 模型 
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具有 更 小 的 残 差 。 
isr " 1.5 b 1.5 : 
1.0 上 1.0 1.0 
0.5 r- 0.5 0.5 
0r 0 E 0 = = 
-0.5 H -0.5 70.5 
-1.0 上 -]1.0 -1.0 
8 
S-15 -]1.5 -]1.5 
$ LR GWLR LR GWLR LR GWLR 
38 1 
5r 1.5 1.5 
R d e f 
1.0 - 1.0 1.0 
05 - 0.5 0.5 
0r n - 0 0 
-0.5 L -0.5 -0.5 
-1.0 - -1.0 -1.0 
-]1.5 -]1.5 -1.5 
LR GWLR LR GWLR LR GWLR 


El4 模型 残 差 箱 形 图 
Fig.4 Box figure the.model residuals 
全 样本 模型 的 残 差 空间 自 相关 (图 5) 结 果 显 示 GWLR 模型 的 Moran's I ff F6 LR 模型 小 ,说 明 与 LR 模型 
相 比 ,GWLR 模型 在 福建 省 林 火 发 生 与 否 的 判别 问题 中 考虑 了 更 多 的 空间 自 相 关 问 题 ,GWLR 模型 对 福建 省 
林 火 发 生 的 拟 合 效果 较 好 。 


气候 条 件 的 变化 对 森林 火灾 的 发 生 具 有 重要 的 影 
响 , 因 此 ,本 文 基 于 气象 因 隆 应 用 传统 的 逻辑 斯 蒂 回 归 
模型 和 地 理 加 权 逻 辑 斯 带 回 归 模 型 对 福建 省 2000 一 
2005 年 的 林 火 数据 进行 分 析 , 对 两 种 模型 的 预测 能 
以 及 对 模型 在 森林 火灾 应 用 方面 的 适用 性 进行 分 析 。 
模型 变量 选择 结果 显示 ,日 最 高 地 表 气 温 "“ 日 最 低 
地 表 气 温 ”“ 目 平均 风速 "“24 小 时 降水 量 ” “日 最 高 
KIVE” H 照 时 数 "“ 日 最 高 气温 ”和 “日 最 低 相 
对 湿度 ”等 8 个 变量 均 是 LR 模型 和 GWLR 模型 的 全 样 0204 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
REPE D S 个 气象 因子 是 福建 省 森林 火灾 的 0 ia CE 
主要 影响 因子 。 和 气温 的 变化 会 使 可 燃 物 的 湿度 发 生变 图 5 全 样本 模型 的 残 差 空间 自 相关 
化 ,改变 火灾 的 气候 条 件 ,从 而 对 森林 火灾 产生 影响 。 Fig.5 Final sample of residual spatial autocorrelation 
我 们 的 研究 结果 显示 “气温 ”对 和 森林 火灾 的 发 生 具 有 
重要 的 影响 ,这 与 Liu 7 A Hu 和 Zhou” 的 研究 结果 一 致 ;相对 湿度 是 反应 林内 可 燃 物 含水 量 的 一 个 重要 指 
标 ,对 林 火 发 生 有 重要 影响 。 我 们 的 研究 结果 表明 “相对 湿度 " 对 林 火 发 生 有 重要 的 影响 ,与 Zhang! 的 研究 
结果 一 致 。 通 常 一 个 地 区 重 特大 火灾 的 发 生 ,与 最 大 风速 是 相关 ,但 是 ,本 文 在 对 气象 因子 进行 分 析 时 ,逻辑 


兰 指数 
Global Moran's I 


全 局 莫 
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斯 带 回归 模型 却 将 “日 最 大 风速 "剔除 ,而 地 理 加 权 人 逻辑 斯 带 回 归 模 型 回归 模型 则 显示 “日 最 大 风速 "是 一 个 
空间 非 平稳 变量 ,因此 ,模型 变量 的 选择 是 否 理 想 与 模型 的 选择 相关 。 此 外 ,本 研究 并 没有 考虑 地 形 植被 类 
型 a 因此 可 能 会 对 研究 结果 产生 偏差 。 
传统 的 逻辑 斯 带 回归 模型 假设 空间 变量 都 是 平稳 变量 ,忽略 了 模型 变量 的 空间 异 质 性 ,模型 的 拟 合 结 

无 法 全 面 反 映 变 量 的 空间 关系 ， Ho m S RET fr EE 3 PK AC EHE RS CAR DSL THU H S iS zs DR] 
ap qa e 与 传统 的 逻辑 斯 带 回归 模型 相 比 ,地 理 加权 逻 辑 斯 蒂 回 归 模 型 的 拟 合 效 果 较 好 这 
与 Saefuddin ^" , Koutsias | Wu 和 Zhang [3 等 前 人 的 研究 结果 一 致 , 表 明 在 对 具有 空间 结构 的 数据 进行 拟 合 
分 析 时 ,应 考虑 空间 地 理 位 置 变化 对 因 变 量 结果 的 影响 。 


4 结论 


本 文 应 用 逻辑 斯 蒂 回 归 模型 和 地 理 加 权 导 和 辑 斯 蒂 回 归 模 型 对 福建 省 森林 火 帘 空间 格局 和 气象 影响 因子 
进行 分 析 ,并 对 两 种 模型 的 拟 合 效 果 进 行 研究 ,结果 表明 :(1) 日 最 高 地 表 气 温 " “日 最 低地 丧气 温 ” “日 平 
均 风 速 "“24 小 时 降水 量 " “日 最 高 本 站 气压 “日照 时 数 ” “日 最 高 气温 ” 和/ 日 最 低 相 对 湿度 "等 8 个 变 
量 是 3 个 模型 的 共同 变量 ,是 影响 福建 省 林 火 发 生 的 主要 气象 因子 ; (2) 在 福建 省 灯火 发 生 的 分 类 判别 中 ,与 
传统 的 逻辑 斯 蒂 回 归 模 型 相 比 ,地 理 加 权 逻 辑 斯 带 回 归 模 型 的 拟 合 效 果 更 好 ,更 适合 福建 省 和 森林 火灾 的 数据 
结构 。 


参考 文献 ( References ) : 


[1] Grogan P, Burns T D, Chapin F S. Fire effects on ecosystem nitrogen cycling in a Californian bishop pine forest. Oecologia, 2000, 122(4): 
537-544. 

[2] Liu H P, Rerson J T, Lindfors J, Chapin F S. Changes in the surface energy budget after fire in boreal ecosystems of interior Alaska; An annual 
perspective. Journal of Geophysical Research Atmospheres, 2005, 110(D13) ; D13101. 

[3] Chas-Amil M L, Touza J, Prestemon P. Spatial distribütion of human-caused forest fires in Galicia (NW Spain). Ecology and the Environment, 
2010, 137; 247-258. 

[4] Bisquert M, Caselles E, Sánchez J M, Casellés. Application of artificial neural networks and logistic regression to the prediction of forest fire danger 
in Galicia using MODIS data. International Journal of Wildland Fire, 2012, 21(8) : 1025-1029. 

[5] Parisien M A, Snetsinger S, Greenbérg J A, Nelson C R, Schoennagel T, Dobrowski S Z, Moritz M A. Spatial variability in wildfire probability 
across the western United States: Intemational Journal of Wildland Fire, 2012, 21(4) ; 313-327. 

[6] Renard Q, Pélissier R, Ramesh B-R,, Kodandapani N. Environmental susceptibility model for predicting forest fire occurrence in the Western Ghats 
of India. International Journal of Wildland Fire, 2012, 21(4) ; 368-379. 

[7] Lozano F J, Suárez-Seoane S, Kelly M, Luis E. A multi-scale approach for modeling fire occurrence probability using satellite data and 
classification trees; A case study in a mountainous Mediterranean region. Remote Sensing of Environment, 2008, 112(3) : 708-719. 

[8] SU, OC, 马 祥 庆 , ARE, 孙 龙 , 胡 海 清 , 杨 婷 婷 . 大 兴安 岭 塔 河 地 区 雷击 火 发 生 驱 动因 子 综合 分 析 . 生态 学 报 , 2015, 35(19) : 
6439- 6448. 
[9] BE, 胡 海 浦 ， 金 森 , 马 志 海 , 张扬 . 基于 负 二 项 和 零 膨 胀 负 二 项 
2010, 34(5) : 571-577. 

[10] RE, WEE, MEM, WE, REE, RE, 潘 建 峰 . 大 兴安 岭 塔 河 地 区 林 火 发 生 的 优势 预测 模型 选择 . 应 用 生态 学 报 , 2014, 25 
(3) : 731-737. 
[u] 郭 福 涛 ,明海 清 , 马 志 海 , 张扬 . 不 同 模型 对 拟 合 大 兴安 岭 林 火 发 生 与 气象 因素 关系 的 适用 性 . 应 用 生态 学 报 , 2010, 21(1) : 159-164. 


[12] Koutsias N, Martínez J, Chuvieco E, Alligówer B. Modeling wildland fire occurrence in southern europe by a geographically weighted regression 


归 模 型 的 大 兴安 岭 地 区 雷击 火 与 气象 因素 的 关系 . 植物 生态 学 报 ， 


n 


approach // Proceedings of the 5th International Workshop on Remote Sensing and GIS Applications to Forest Fire Management; Fire Effects 
Assessment. Spain; Universidad de Zaragoza, 2005; 57-60. 

[13]  Martínez-Fernández J, Chuvieco E, Koutsias N. Modelling long-term fire occurrence factors in Spain by accounting for local variations. with 
geographically weighted regression. Natural Hazards and Earth System Sciences, 2013, 13(2) : 311-327. 


[14] Rodrigues M, De La Riva J, Fotheringham S. Modeling the spatial variation of the explanatory factors of human-caused wildfires in Spain using 


http ://www.ecologica.cn 


201706.00841v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


4140 E dm 学 报 37 卷 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[33] 


geographically weighted logistic regression. Applied Geography, 2014, 48; 52-63. 

Wang Q, Ni J, Tenhunen J. Application of a geographically-weighted regression analysis to estimate net primary production of Chinese forest 
ecosystems. Global Ecology and Biogeography, 2005, 14(4) : 379-393. 

Tu J. Spatially varying relationships between land use and water quality across an urbanization gradient explored by geographically weighted 
regression. Applied Geography, 2011, 31(1) ; 376-392. 

Propastin P. Modifying geographically weighted regression for estimating aboveground biomass in tropical rainforests by multispectral remote sensing 
data. International Journal of Applied Earth Observation and Geoinformation, 2012, 18; 82-90. 

Freedman R, Roy S S. Spatial patterning of Manta birostris in United States east coast offshore habitat. Applied Geography, 2012,-32(2); 
652-659. 

Aguilar G D, Farnworth M J. Stray cats in Auckland, New Zealand; discovering geographic information for exploratory spatial analysis. Applied 
Geography, 2012, 34; 230-238. 

国家 林业 局 . 第 八 次 全 国 森 林 资 源 清查 结果 . 林业 资源 管理 , 2014, (1) : 1-2. 

李 兆 明 , AJEET. 福建 森林 火险 天 气 等 级 预测 预报 方法 . 福建 林学 院 学 报 , 1989, 9(2) : 172-176. 

朱 学 平 , 何 东 进 , 丁 福 立 , 洪 伟 , WOR, 纪 志 荣 . 基于 突变 级 数 法 的 福建 省 森林 火灾 评价 分 析 . 福建 林学 院 学 报 ，2011, 31(4) : 
295-299. 


Garcia C V, Woodard P M, Titus S J, Adamowicz W L, Lee B S. A logit model for predicting the daily occurrence. of human caused forest-fires. 
International Journal of Wildland Fire, 1995, 5(2) ; 101-111. 

Zhang Z X, Zhang H Y, Zhou D W. Using GIS spatial analysis and logistic regression to predict the probabilities of human-caused grassland fires. 
Journal of Arid Environments, 2010, 74(3) : 386-393. 

del Hoyo L V, Isabel M P M, Vega F J M. Logistic regression models for human-caüsed wildfire Tisk estimation; analysing the effect of the spatial 
accuracy in fire occurrence data. European Journal of Forest Research, 2011 ,.130(6).; 983-996. 

郭 含 站 , SK EE, 徐 丽华 , RIRE, KAK. 基于 地 理 加 权 回 归 的 区 域 森林 碳 储量 估计 . 浙江 农林 大 学 学 报 , 2015, 32(4) : 497-508. 
HRG, 田 铮 , AEs. 时 空地 理 加 权 回 归 模 型 的 时 空 非 平稳 性 检验 ”系统 国 程 理论 与 实践 , 2013, 33(6) : 1537-1542. 

Wu W, Zhang L J. Comparison of spatial and non-spatial logistic regression models for modeling the occurrence of cloud cover in northeastern Puerto 


Rico. Applied Geography, 2013, 37; 52-62. 


n 


Zhang H J, Qi P C, Guo G M. Improvement of fire dangér modelling with geographically weighted logistic model. International Journal of Wildland 
Fire, 2014, 23(8) ; 1130-1146. 

Rodrigues M, De La Riva J. An insight into machine-learning algorithms to model human-caused wildfire occurrence. Environmental Modelling & 
Software, 2014, 57; 192-201. 

Saefuddin A, Setiabudi N A, Fitrianto A. On. comparison between logistic regression and geographically weighted logistic regression; with 
application to Indonesian poverty data. World Applied Sciences Journal, 2012, 19(2) ; 205-210. 

Liu Z H, Yang J, Chang Y „Weisberg P J, He H S. Spatial patterns and drivers of fire occurrence and its future trend under climate change in a 
boreal forest of Northeast China. Global Change Biology, 2012, 18(6) : 2041-2056. 

Hu T Y, Zhou G S. Drivers of lightning-and human-caused fire regimes in the Great Xing'an Mountains. Forest Ecology and Management, 2014, 
329; 49-58. 


http ://www.ecologica.cn 


ChinaXiv 合 作 期 刊 


12 期 梁 慧 玲 ”等 :比较 逻辑 斯 蒂 与 地 理 加 权 逮 辑 斯 蒂 回归 模型 在 福建 林 火 发 生 的 适用 性 4141 
MR: 
附 表 1 LR 模型 和 GWLR 模型 的 系数 估计 (样本 1) 
Annexed table 1 LR model and GWLR model coefficient estimation ( sample 1) 
"Tm — 日 最 高 地 日 最 低地 日 平均 最 大 24 小 时 日 最 高 本 ec 日 平均 日 最 小 
型 统计 量 BEJ i à ^ M 5 . H KEI 日 最 高 气温 
i e BETO dedi Rb a MO Rok Gg PANA RANTE antep — quip 
i xd n GST. max GST. min Wnd avg Wnd, max Pre Pres max dm RH avg RH min 
LR “估计 系数 -3.979183 0.004932 -0.020798 0.011638 -0.000071 -0.004711 0.000243 0.012919 0.007426 -0.003009 -0.003736 
标准 误差 1.064268 0.000661 0.001352 0.004287 0.002749 .001232 0.000096 0.001591 0.001519 0.005586 0.005317 
个 计 系 数 -标准 误差 -5.043451 — 0.004271 -0.02215 0.007351 -0.00282 -0.005943 — 0.000147 0.011328 0.005907 -0.008595  =0:009053 
估计 系数 + 标准 误差 -2.914915 0.005593 -0.019446 0.015925 0.002678 -0.003479 0.000339 0.014510 0.008945 0.002577 0.001581 
GWLR 最 小 值 -54.198396 -0.002282 -0.033631 -0.065206 -0.022832 -0.017348 -0.000429 -0.001035 -0.000664 =0.0704220 ì} -0.037178 
四 分 之 一 分 位 数 -11.951030 0.002654 -0.026593 -0.000398 -0.008862 -0.006708 0.000091 0.008001 0.006217 -0.018073 -0.013137 
平均 数 -9.182290 0.003929 -0.023161 0.008436 0.001044 -0.005226 — 0.000746 0.011935 0.011610 -0.003814 -0.004102 
中 位 数 -6.484965 0.003700 -0.023076 0.012048 0.001753 -0.004936 0.000527 0.013266 0.010479 0.003381 -0.005216 
= 四 分 之 三 分 位 数 -2.419304 0.004719 -0.019046 0.024120 0.010857 -0.003099 0.000921 0.015655 ~ | 0.015613 0.013657 0.003580 
p 
= 最 大 值 3.108987 — 0.010752 -0.015394 0.061361 0.027887 .000520 0.005137 — 0.022256 40.031157 0.029408 0.037951 
T- 四 分 位 数 范围 9.531726 0.002065 0.007547 0.024518 0.019719 .003609 A 0.00083 0.007654 — 0.009396 0.08173 0.016717 
"T 
cO 
= 
«, 
Te 附 表 2 GWLR 模型 空间 ( 非 ) 平稳 变量 检验 结果 
hel 
= Annexed table 2 GWLR model of space ( nor) smooth variable inspection results 
N 变量 样本 1 样本 2 样本 3 样本 4 样本 5 最 终 样 本 样本 显著 个 数 
YT Variable Sample 1 Sample2 — Sample 3  Sample4 — Sample 5 Final sample No.samples signif 
pm 
© 日 最 高 地 表 气 温 
N . 十 T 十 十 十 十 5 
- The highest surface temperature 
Lu E E 
9 2 最 低 电表 气温 + 十 十 十 十 十 5 
Su The minimum surface temperature 
>< 日 平均 风速 
a . : 十 十 十 十 十 十 5 
x Daily average wind speed 
[( vE 
B 日 最 大 风速 
em Dail i nnd speed " * * * + + 5 
c aily maxium wind spee 
© 24 小 时 降水 量 + + + + + + 3 
Daily precipitation 
日 最 高 本 站 气压 
. . 十 十 十 十 十 十 5 
The highest station pressure 
Ha mp 
日 照 时 数 十 十 十 十 十 十 5 
Hours of sunshine 
E Pr yH 
H TIR A i 十 十 十 十 十 十 5 
Daily maxunum temperature 
日 平均 相对 湿度 
十 十 十 十 十 十 5 
Daily average relative humidity 
日 最 小 相对 湿度 
十 十 一 十 十 十 4 


The minimum relative humidity 


4 表示 该 变量 是 空间 非 平稳 变量 ;- 表 示 该 变量 是 空 上 
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